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We perform a global analysis of neutrino oscillation data, including high-precision measurements 
of the neutrino mixing angle #13 at reactor experiments, which have confirmed previous indications 
£SJ ' in favor of #13 > 0. Recent data presented at the Neutrino 2012 Conference are also included. We 

t-H ' focus on the correlations between #13 and the mixing angle 623, as well as between 613 and the 

neutrino CP- violation phase S. We find interesting indications for 623 < 7r/4 and possible hints for 
CNJ ■ S ~ 7r, with no significant difference between normal and inverted mass hierarchy. 

c : 
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I— 5 ' I. INTRODUCTION 

00 ■ 

Current neutrino oscillation experiments (except for a few anomalous results) can be interpreted in a simple three- 
neutrino framework, where the three flavor states v a — {v ei u^,v T ) are quantum superpositions of three light mass 
states Vi = (yi, v 3 , z^) via a unitary mixing matrix U a i, parametrized in terms of three mixing angles (#12, 013, $23) 
O !■ and one possible CP-violating phase 5 in standard notation 0, Q ■ 

In neutrino oscillations, CP violation is a genuine 3v effect which may be observed (provided that S 7^ 0, ir) only if 
', all the mixings 0,j and the squared mass differences m 2 — m 2 are nonzero @ . The latter condition is experimentally 
established, and can be expressed in terms of the two independent parameters Sm 2 = m 2 — ra\ > 1] and Am 2 = 
l— ~ '] m 2 — (m 2 + m 2 )/2 [|[, where Am 2 > (< 0) corresponds to normal (inverted) mass spectrum hierarchy. 

Until very recently, the further condition Oij ^ could be considered as established for 612 and #23 [J , and quite likely 
(at ~ 3(7 level) but not conclusively settled for #13 This year, the short-baseline (SBL) reactor experiments Daya 
Tj" ■ Bay Q and RENO [?J have definitely established that 6*13 > at ~ 5a, by observing T7 e disappearance from near to far 
IT) ; detectors. In particular, Daya Bay and RENO have measured sin 2 6 13 ~ 0.023 ± 0.003 @ and sin 2 6 13 ~ 0.029 ± 0.006 
O^l 1 @; 1 respectively. Consistent indications were also found in the Double Chooz reactor experiment with far detector 
only (sin 2 #13 ~ 0.028 ± 0.010) [l(l 03 ■ All these reactor data are in good agreement with the results of our latest 
I/"") I global analysis of oscillation data in 5], which provided sin 2 6*13 = 0.021-0.025 at best fit, with a la error of ±0.007. 

It should be remarked that we had previously obtained hints in favor of sin 2 13 ~ 0.02 from a detailed analysis of 
solar and long-baseline reactor data [l2|,[l3| (see also [l4{ for similar, independent hints), consistently with an earlier 
(weak) preference for 6*13 > from atmospheric neutrinos [1, [l3[ . The hints became a ~2cr indication for 6*13 > in 
combination with early appearance data from the MINOS long-baseline accelerator experiment fl5| . and provided a 
> 3a evidence by including the remarkable low-background appearance data from the T2K experiment [j| . The Daya 
Bay and RENO measurements have shown that our global 3v analyses in (g l lij 13| — the latest of a series started 
two decades ago [lj| — were on the right track in the hunt to #13. See also 17 for other recent analyses of #13 
constraints prior to the Daya Bay and RENO results. 

With sin 2 #13 as large as 2-3 xlO -2 , the door is open to CP violation searches in the neutrino sector, although the 
road ahead appears to be long and difficult [2(| HH- At present, it makes sense to squeeze, from the available data, 
any tiny bit of information about S. An interesting attempt has been made in [22| , using reactor and accelerator data. 
However, atmospheric v data may also usefully probe S 0, HH. To this purpose, we update the analysis in [5] by 
including new atmospheric, LBL accelerator and SBL reactor data, as available after the Neutrino 2012 Conference 
i- We have also extended our atmospheric v codes (previously limited to cos S = ±1 0,@) to generic values of 8. 
Among the results obtained, we pay particular attention on a possible preference in favor of #23 < 7r/4 and of 5 ~ tt 
in both hierarchies (although with limited statistical significance). We also discuss the implications of the oscillation 
parameter constraints for absolute v mass searches, as well as some limitations and challenges of global analyses. 

The present work is structured as follows. In Sec. II we describe some methodological issues, which may be skipped 
by readers interested only in the main results. In Sec. Ill we discuss the results of our analysis in terms of covariance 
among the parameters (sin 2 #13, sin 2 #237 S), for both normal and inverted hierarchy. In Sec. IV we summarize the 
constraints on the mass-mixing oscillation parameters, and describe their implications for the observables sensitive 
to absolute neutrino masses. We conclude our work in Sec. V. Details of atmospheric neutrino flavor evolution for 
generic 5 are confined in the Appendix. 
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II. METHODOLOGY: GROUPING AND ANALYZING DIFFERENT DATA SETS 

No single oscillation experiment can sensitively probe, at present, the full parameter space spanned by 
(5m 2 , ±A?n 2 , 612, $13, $23, 5). Therefore, it is necessary to group in some way the experimental data, in order 
to study their impact on the oscillation parameters. For instance, in |5J we showed that consistent indications in favor 
of nonzero #13 emerged from two different datasets, one mainly sensitive to 5m 2 (solar plus KamLAND experiments) 
and another mainly sensitive to Am 2 (CHOOZ plus atmospheric and LBL accelerator experiments). In this work 
we adopt an alternative grouping of datasets, which is more appropriate to discuss interesting features of the current 
data analysis, such as the covariance among the parameters (sin 2 #13, sin 2 823, 5) in both mass hierarchies. 

A. LBL + solar + KamLAND data 

We remind that LBL accelerator data (from the K2K, T2K, and MINOS experiments) in the — > disappearance 
channel probe dominantly the Am 2 -driven amplitude 

l^3| 2 (l" l^al 2 ) =COS 2 1 3Sin 2 23 (l-COS 2 1 3Sin 2 2 3) , (1) 

which is slightly octant-asymmetric in 623 for #13 ^ 0. In the — > v e appearance channel, the dominant Am 2 -driven 
amplitude is 

|t^ 3 | 2 |t4 3 | 2 = cos 2 9 13 sin 2 9 13 sin 2 9 23 , (2) 

which is definitely octant-asymmetric in 623 for 813 7^ 0. In both the appearance and the disappearance channels, 
subdominant terms driven by Sm 2 and by matter effects can also contribute to lift the octant symmetry and to 
provide some weak sensitivity to sign(Am 2 ) and to 5, see e.g. [24j for a recent perturbative approach at "large" 6*13. 
As already noted in Q , the T2K and MINOS indications in favor of — » v e appearance induce an anti-correlation, 
via Eq. between the preferred values of sin 2 623 and sin 2 $13. This covariance is relevant in the analysis of the 623 
octant degeneracy [25| and has an indirect impact also on the preferred ranges of 5 via subdominant effects. 

In order to make the best use of LBL accelerator data, it is thus useful to: (1) analyze both disappearance 
and appearance data at the same time and in a full 3v approach; (2) combine LBL with solar and KamLAND 
data, which provide independent constraints on (5m 2 , 812, #13) and thus on the subdominant 7>v oscillation terms. As 
discussed below, once the (relatively well known) oscillation parameters sin 2 6*12, 5m 2 and Am 2 are marginalized away, 
interesting correlations emerge among the remaining parameters (sin 2 #13, sin 2 $23, 5). Conversely, these interesting 
bits of information are partly lost if LBL disappearance data are analyzed in the 2v approximation and/or separately 
from appearance data, as it has often been the case in official analyses by experimental collaborations. 

In this work, the previous LBL data used in Q are updated with the inclusion of the first T2K disappearance 
constraints [2(| and of the latest T2K appearance data [27[. We note that recent MINOS disappearance data [28| 
are no longer in disagreement with previous results. Therefore, it makes sense to use both v and V MINOS disap- 
pearance constraints, which we take from [29| . together with updated MINOS appearance data. For later purposes, we 
note that recent T2K and (especially) MINOS data are best fit for slightly nonmaximal mixing (sin 2 26> 23 ^ 0.94-0.98 
[26l [28l [29j ) roughly corresponding to the octant-symmetric values sin 2 #23 ~ 0.4 or 0.6). A slight preference for 
nonmaximal mixing emerged also from our analysis of K2K LBL data in Q . 

B. Adding SBL reactor data 

After grouping LBL accelerator plus solar plus KamLAND data (LBL + solar + KamLAND), it is important to add 
the independent and "clean" constraints on #13 coming from SBL reactor experiments in the v e — ¥ v e disappearance 
channel, which probe dominantly the Am 2 -driven amplitude 

\Ue3?(l ~ \Ue3?) = sin 2 6 13 cos 2 6 13 . (3) 

In the reactor dataset, subdominant terms are slightly sensitive to (5m 2 , 612) and, as noted in [30| and discussed in 
[3l| . probe also the neutrino mass hierarchy. We include far-detector data from CHOOZ |32| and Double Chooz [TTJ] 
and near-to-far detector constraints from Daya Bay Q and RENO 0, E|. We do not include data from pre-CHOOZ 
reactor experiments, which mainly affect normalization issues. 

Indeed, the analysis of reactor experiments without near detectors depends, to some extent, on the absolute nor- 
malization of the neutrino fluxes, which we choose to be the "old" (or "low") one, in the terminology of j5|. We shall 
also comment on the effect of adopting the "new" (or "high" ) normalization recently proposed in [33, [34| . Constraints 
from Daya Bay and RENO are basically independent of such normalization, which is left free in the official analyses 
and is largely canceled by comparing near and far rates of events 0, 0|- At present, it is not possible to reproduce, 
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from published information, the official Daya Bay and RENO data analyses with the permill accuracy appropriate to 
deal with the small systematics affecting near/far ratios. We think that, for the purposes of this work, it is sufficient 
to take their measurements of sin 2 2#i3 at face value, as gaussian constraints on such parameter. Luckily, such con- 
straints appear to depend very little on the Am 2 parameter within its currently allowed range; see the (Am 2 , sin 2 26*13) 
prospective sensitivity plots in [35| (Daya Bay) and [36| (RENO). Of course, a joint analysis of all SBL reactor data 
made by the current collaborations would be desirable, since a few systematics are correlated among the experiments. 



As shown in [25|, LBL data in disappearance and appearance mode generally select [via Eqs. ([T]) and ©], two 



degenerate (#23, $13) solutions, characterized by nearly octant-symmetric values of #23 and by slightly different values 
of #13. By selecting a narrow range of $13, precise reactor data can thus (partly) lift the #23 octant degeneracy J2a 
(see also [37|)- Amusingly, the fit results in Sec. Ill resemble the hypothetical, qualitative iv scenario studied in [25j . 



C. Atmospheric neutrino data 



After combining the (LBL + solar + KamLAND) and (SBL reactor) datasets, we finally add the Super-Kamiokande 
atmospheric neutrino data (SK atm.), as reported for the joint SK phases I-IV in (38| (but with no statistical vjv 
separation [38| , which we cannot reproduce in detail) . The SK data span several decades in neutrino and antineutrino 
energy and pathlengths, both in vacuum and in matter, in all appearance and disappearance channels involving 
and v e , and thus they embed an extremely rich 3is oscillation physics. 

In practice, it is difficult to infer — from atmospheric data — clean Zv information beyond the dominant parameters 
(Am 2 , 023 )• Subdominant oscillation effects are often smeared out over wide energy-angle spectra of events, and can 
be partly mimicked by systematic effects. For this reason, "hints" coming from current atmospheric data should 
be taken with a grain of salt, and should be possibly supported by independent datasets. For instance, we have 
attributed some importance to a weak preference for #13 > found from atmospheric SK data in Q , only after it was 
independently supported by solar+KamLAND data [13| and, later, by LBL accelerator data Q. Similarly, we have 
typically found a preference of atmospheric SK data for #23 < 7r/4 [J, Q; in the next Section, we shall argue that such 
preference now finds some extra support in other datasets, and thus starts to be an interesting frontier to be explored. 

The situation is more vague for 8. We argued in [4| (and also found in [5|) that a slight electron excess in the 
atmospheric event spectra at sub-GeV energies could be better fit with cos (5 = — I as compared with cos 8 = +1, 
via interference terms [3, [23| in the oscillation probability. Since the analyses in 0, 0] were limited to the two 
CP-conserving cases cos 8 = ±1, we have now extended our atmospheric neutrino codes to generic values of 8 in the 
oscillation probability; details are given in the Appendix. We continue to find a preference for cos 8 ~ —1, as described 
in the next Section. This possible hint for 8 ~ ir is roughly consistent with the SK official (although preliminary) 
analyses in [38|, 39] , but is not clearly matched by a similar hint coming from other data. This is another reason for 
choosing to present atmospheric constraints only after the discussion of other datasets. In conclusion, we think that 
is methodologically useful to show, in sequence, the impact of data from (LBL + solar + KamLAND), plus (SBL 
reactors), plus (SK atm.) experiments on the neutrino oscillation parameters. 



D. Limitations and challenges of global analyses 



Our global analyses offer contributions to the discussion on the neutrino oscillation phenomenology, but should not 
be considered as a substitute for the official oscillation analyses performed by the experimental collaborations, which 
include unpublished or unreproduciblc information. Therefore, our estimated parameter ranges may be slightly offset 
with respect to those estimated by the collaborations in dedicated 3v data analyses (when available). Our educated 
guess is that possible offsets are < la at present, and often much lower. Of course, even a fraction of one standard 
deviation may matter when discussing hints at or below the 2a level, as done in the next Section. However, the 
success story of the indications of #13 > [!, [l3| shows that discussions of ~ 2a effects may still have some interest. 

Global 3v analyses will face several new challenges in the near future. As already remarked, a joint analysis 
of all reactor data with near and far detectors (Daya Bay, RENO, Double Chooz) will be useful to get the most 
stringent constraints on $13. The T2K and MINOS long-baseline accelerator experiments are urged to abandon any 
2v approximation in the interpretation of their (disappearance) data, and focus on full-fledged 3is combinations of 
appearance plus disappearance data. Increasing attention should be paid to refined features of the LBL analysis, such 
as the impact of cross section assumptions on the oscillation parameter ranges [40]. Future solar and long-baseline 
reactor data might slightly reduce the uncertainties of the (#12, <5m 2 ) parameters, which drive subleading oscillation 
terms at higher energies. Concerning atmospheric ^s and their associated systematics, we think that, while waiting 
for future large-volume detectors and data, the existing SK atmospheric data have not yet exhausted their physics 
potential: dedicated 2>v analyses from the SK collaboration might reveal intriguing indications on #23 and on 8, 
especially if their Monte Carlo simulations were reprocessed by assuming full, unaveraged 3z/ oscillations from the 
very beginning (rather than re- weighting unoscillated simulations with factors embedding averaged oscillations 41]). 
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III. RESULTS: CORRELATIONS BETWEEN 6»i 3 , 23 AND 5 

In this section we focus on two emerging features of our analysis: converging hints in favor of #23 < 7r/4, and a 
possible (weak) hint in favor of 5 ~ it. The correlations of 823 and <5 with #13 are discussed in some detail. As in our 
previous works allowed regions are shown at No confidence levels, where No = y/ A\ 2 [lj]. It is understood 

that, in each figure, undisplayed oscillation parameters have been marginalized away. 

Figure 1 shows the results of the analysis in the plane (sin 2 6*13, sin 2 #23), for both normal hierarchy (NH, upper 
panels) and inverted hierarchy (IH, lower panels). From left to right, the panels refer to increasingly rich datasets: 
LBL accelerator + solar + KamLAND data (left), plus SBL reactor data (middle), plus SK atmospheric data (right). 

In the left panels, LBL appearance data anti-correlate sin 2 #13 and sin 2 623 via Eq. ([2]). On the other hand, LBL 
disappearance data (via their current preference for sin 2 2#23 < 1) disfavor maximal mixing at > lo. As a consequence, 
two quasi-degenerate \ 2 minima emerge at complementary values of sin 2 #23 and at somewhat different values of #13. 
The degeneracy is slightly lifted by solar+KamLAND data, whose preference for sin 2 #13 ~ 0.02 [5fl picks up the first 
octant solution in NH, and the second octant solution in IH. However, as far as LBL+solar+KamLAND data are 
concerned, the statistical difference between the two #23 solutions remains negligible (< 0.3<r) in both NH and IH. 

In the middle panels, the addition of SBL reactor data (most notably from Daya Bay and RENO) fixes sin 2 #13 
with high accuracy and at relatively "large" values, which are best matched at low #23 — hence the overall preference 
for the first #23 octant in both hierarchies. Such preference is more pronounced in NH (at the level of ~ lo). In IH, 
both T2K and MINOS appearance data can accommodate values of 813 generally larger than in NH [53, H3, IH, HH 
(as also evident from the left panels), so that the agreement with SBL reactor data can be easily reached in both 
octants, with only a small preference (~0.4cr) for the first. The combination of LBL accelerator and SBL reactor data 
to lift the octant degeneracy was proposed in |25[ . 

In the right panels, atmospheric v data do not noticeably improve the constraints on #13, but corroborate the 
preference for the first octant (as already found in [HQ), in both NH (slightly below the 3o level) and IH (slightly 



below the 2o level). [We do not observe an octant flip with the hierarchy as in [38[.] In conclusion, from Fig. 1 
we derive that both atmospheric and non-atmospheric v data seem to prefer, independently, the first octant of #23 
(especially in normal hierarchy), with a combined statistical significance < 3<r in NH and < 2o in IH. 




FIG. 1: Results of the analysis in the plane charted by (sin 2 #13, sin 2 #23), all other parameters being marginalized away. From 
left to right, the regions allowed at 1, 2 and 3a refer to increasingly rich datasets: LBL+solar+KamLAND data (left panels), 
plus SBL reactor data (middle panels), plus SK atmospheric data (right panels). Best fits are marked by dots. A preference 
emerges for 6*23 in the first octant in both normal hierarchy (NH, upper panels) and inverted hierarchy (IH, lower panels). 
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Figure 2 shows the results of the analysis in the plane (sin 2 6*13, 8). The conventions used are the same as in Fig. 1. 
Since the boundary values 8/tt = and 2 are physically equivalent, each panel could be ideally "curled" by smoothly 
joining the upper and lower boundaries. 

In the left panels, constraints on sin 2 #13 are placed both by solar+KamLAND data (independently of 8) and by 
current LBL accelerator data (somewhat sensitive to 8). Once more, it can be noted that larger values of 6*13 are 
allowed in IH. The best fit points are not statistically relevant, since all values of 8 provide almost equally good fits 
at ~ la level. The "fuzziness" of the la contours is a consequence of the statistical degeneracy of the two solutions 
allowed at la in Fig. 1, and which involve complementary values of 823 and somewhat different values of #13. At la, 
the fit is "undecided" between the wavy bands at smaller and larger values of 6*13 , and easily flips between them. At 
2 or 3a the two bands merge and such degeneracy effects are no longer apparent. 

In the middle panels, SBL reactor data pick up a very narrow range of 13 and suppress degeneracy effects. Some 
sensitivity to 8 starts to emerge, since the "wiggles" of the bands in the left panel best match the (5-independent SBL 
reactor constraints on sin 2 813 only in certain ranges of 8. The match is generally easier in inverted hierarchy (where 
LBL data allow a larger 6*13 range) than normal hierarchy. 

In the right panels, atmospheric neutrino data induce a preference for 8 ~ 7r, although all values of 8 are still allowed 
at ~ la. Such a preference is consistent with our previous analyses limited to cosi5 = ±1 [J, |5[, where we found 8 = tt 
preferred over 8 = 0, in both normal and inverted hierarchy. As discussed in for 8 ~ tt the interference term in 
the oscillation probability provide some extra electron appearance in the sub-GeV atmospheric neutrino data, which 
helps fitting the slight excess of electron-like events in this sample. In our opinion, atmospheric data can provide 
valuable indications about the phase S, which may warrant dedicated analyses by the SK experimental collaboration, 
especially in combination with data from the T2K collaboration, which uses SK as far detector and thus shares some 
systematics related to final state reconstruction and analysis. 

Concerning the hierarchy, in the middle panels of Figs. 1 and 2 (all data but SK atm.) we find a slight preference 
for IH with respect to NH (A% 2 ~ —0.38). The situation is reversed in the right panels (all data, including SK atm.), 
where NH is slightly favored (A% 2 — +0.35). These fluctuations between NH and IH fits are statistically irrelevant. 
We conclude that, in our analysis of oscillation data, there are converging hints in favor of 6*23 < tt/A (especially in 
NH), a possible hint in favor of 8 ~ tt (mainly from SK atm. data), and no hint about the mass hierarchy. 




FIG. 2: Results of the analysis in the plane charted by (sin 2 613, 8), all other parameters being marginalized away. From left 
to right, the regions allowed at 1, 2 and 3cr refer to increasingly rich datasets: LBL+solar+KamLAND data (left panels), plus 
SBL reactor data (middle panels), plus SK atmospheric data (right panels). A preference emerges for 5 values around tt in 
both normal hierarchy (NH, upper panels) and inverted hierarchy (IH, lower panels). 
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IV. SUMMARY OF OSCILLATION CONSTRAINTS AND IMPLICATIONS FOR ABSOLUTE MASSES 



In this section we summarize the previous results in terms of one-parameter constraints, all the others being 
marginalized away. We also show updated oscillation constraints on the main absolute mass observables l45| . 
namely, the effective electron neutrino mass mp (probed in /3 decay), the effective Majorana mass (probed in 0^2/3 
decay searches), and the sum of neutrino masses E, which can be probed by precision cosmology. 

Figure 3 shows the Na bounds on the 3v oscillation parameters. Blue (solid) and red (dashed) curves refer to 
NH and IH, respectively. The curves are expected to be linear and symmetric around the best fit only for gaussian 
uncertainties. This is nearly the case for the squared mass differences Sm 2 and Am 2 , and for the mixing parameters 
sin 2 6*i2 and sin 2 #13. The bounds on sin 2 #23 are rather skewed towards the first octant, which is preferred at < 2a in 
NH and < 3d in IH. Also the probability distribution of d is highly nongaussian, with some preference for 5 close to 
7r, but no constraint above ~2cr. As expected, there are no visible differences between the NH and IH curves for the 
parameters 5m 2 and sin 2 #12, and only minor variations for the the parameters Am 2 and sin 2 #13. More pronounced 
(but < la) differences between NH and IH curves can be seen for sin 2 #23 and, to some extent, for S. 



Synopsis of global 3v oscillation analysis 





■ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




q 1 1 1 1 1 1 1 1 1 1 yi 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 "a i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
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FIG. 3: Results of the global analysis in terms of Na bounds on the six parameters governing 3v oscillations. Blue (solid) and 
red (dashed) curves refer to NH and IH, respectively. 
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TABLE I: Results of the global 3u oscillation analysis, in terms of best-fit values and allowed 1, 2 and 3<r ranges for the 3v 
mass-mixing parameters. We remind that Am 2 is defined herein as m\ — (ml + m\)/2, with +Am 2 for NH and — Am 2 for IH. 
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Table I reports the bounds shown in Fig. 3 in numerical form. Except for 5, the oscillation parameters are constrained 
with significant accuracy. If we define the average la fractional accuracy as l/6th of the ±3ct variations around the 
best fit, then the parameters are globally determined with the following relative precision (in percent): Sm 2 (2.6%), 
Am 2 (3.0%), sin 2 6l 12 (5.4%), sin 2 13 (10%), and sin 2 9 23 (14%). 

A final remark is in order. As noted in Sec. II B, two alternative choices were used in [5| for the absolute reactor flux 
normalization, named as "old" and "new," the latter being motivated by revised flux calculations. Constraints were 
shown in [5j for both old and new normalization, resulting in somewhat different values of 9\2 and #13. The precise 
near/far data ratio constraints from Daya Bay [a, HI and RENO 0, Q are largely independent of such normalization 
issues, which persists only for the reactor data without near detector (i.e., KamLAND, CHOOZ and Double Chooz 
data in this work), with very small effects on the global fit. For the sake of precision, we remark that the values 
in Table I refer to our fit using the "old" normalization for KamLAND, CHOOZ and Double Chooz. By using the 
"new" normalization, the only noticeable effects would be the following overall shifts, with respect to the numbers in 
Table I: A sin 2 fl^/lO -1 ^ +0.05 and A sin 2 6>i 3 /10~ 2 ~ +0.08 (i.e., at the level of ~ 1/3 of a standard deviation). 




£ (eV) m,, (eV) 



FIG. 4: Constraints induced by oscillation data (at 2a level) in the planes charted by any two among the absolute mass 
observables m$ (effective electron neutrino mass), mpp (effective Majorana mass), and E (sum of neutrino masses). Blue (red) 
bands refer to normal (inverted) hierarchy. 
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Let us now discuss the interplay of oscillation and nonoscillation data. The constraints in Table Q] induce strong 
covariances among the three main observables which are sensitive to the absolute masses, namely, mp, mpp and S 
(see EH for notation). Figure 4 shows such covariances in terms of 2a constraints (bands) in the planes charted 
by any couple of the absolute mass observables. As compared to previous results [13, EH, the bands in the {nip, S) 
plane of Fig. 4 are narrower, due to the higher accuracy reached in the determination of all the oscillation parameters. 
Note that, in principle, precise measurements of {mp, S) in the sub-eV range (where the bands for NH and IH branch 
out) could determine the mass spectrum hierarchy. In the two lower panels of Fig. 4, there remains a large vertical 
spread in the allowed slanted bands, as a result of the unknown Majorana phases in the mpp components, which 
may interfere either constructively (upper part of each band) or destructively (lower part of each band). In principle, 
precise data in either the (mpp, mp) plane or the {mpp, E) plane might thus provide constraints on the Majorana 
phases. 

Progress in constraining the neutrino mass and mixing parameters will hopefully lead to a deeper understandin g of 
their origin. Theoretical options range from "accidental" parameter values with no special significance or structure [46[ 
to "special" values pointing towards underlying symmetries [13], just to name a few possibilities in the vast literature 
on models. Precision measurements of neutrinos masses, mixings and phases will provide valuable information to 
narrow this wide theoretical spectrum. 

V. CONCLUSIONS 

We have performed a global analysis of neutrino oscillation data, including recent, high-precision measurements of 
the neutrino mixing angle #13 at reactor experiments (which have confirmed previous indications in favor of #13 > 
[E EH) and updated data released at the Neutrino 2012 Conference @. We have explored the current correlations 
between the mixing parameters sin 2 #13 and sin 2 $23, as well as between sin 2 #13 and the CP-violation phase 5. We 
have found some interesting indications in favor of 6*23 < 7r/4 (at < 3<r in NH and < 2a in IH), as well as possible 
hints of S ~ 7r, but no significant difference between normal and inverted mass hierarchy. We surmise that full-fledged 
2>v analyses of LBL and atmospheric neutrino data by the experimental collaborations would be very useful to better 
assess the statistical relevance of these possible hints. 

Note added. After this work was basically completed, we noted the results of another analysis including recent 
reactor data [Hj]. Some differences with our results emerge in the favored ranges for 623 and 5; they might be due, in 
part, to a different approach to atmospheric neutrino oscillations (which, in our case, do include 5m 2 and 5 effects). 
We also noted the preliminary results of the full 2>v global analysis in [4g|, where #23 < 7r/4 is also preferred. 
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Appendix: Atmospheric neutrino flavor evolution for generic 8 

Atmospheric neutrinos traverse the atmosphere and several Earth shells before being detected. We adopt a five-shell 
approximation of the electron density N in the Earth, in which each j-th cell has sharp edge discontinuities and a 
mild dependence Nj(r) in terms of the normalized radial distance r from the Earth center [50], that can be well 
approximated by a quartic polinomial (5lj 

iV,(r)= aj +/3,r 2 + 7j r 4 , (A.l) 

where the coefficients ay , f3j and jj are given in Table I of |5l[ . 

The evolution operator for atmospheric neutrinos can be written as the product of the evolution operator in each 
shell chord 



7iarth = r (FbK) ■ r (f\p7) 



T {Pm-\Pm) • 7~v (PmPa) , 



(A.2) 



9 



where Pq is the detection point, M the number of shells crossed by neutrinos and Pa the production point in 
atmosphere. The last operator embeds the propagation in atmosphere, governed by the "vacuum" Hamiltonian T~i v . 
Notice that for a real Hamiltonian the calculation of 7Earth can be further simplified using the symmetry properties 
of the electron density along the neutrino path inside the Earth (see appendix B of [5l|). This property is no longer 
valid when the neutrino mixing matrix is not real, i.e., Sep 0, n. 

A first-order approximation for the evolution operator inside the fc-th shell is to consider the electron density 
constant, and equal to the average along the shell chord 

T (iViP fe ) - exp [-i(n v + V k ) ■ D k ] , (A.3) 

where V k = diag{v / 2G , F-/Vfc, 0, 0} is the matter potential, D k the distance travelled by the neutrino inside the shell, 
and 



1 



N k = — I dxN k {^ X 2+ sin 2 rj ) , (A.4) 



where 

r 2 = x 2 + sin 2 77 , (A. 5) 

77 being the nadir angle of the neutrino direction. Handy subroutines for calculating exponentials of real or complex 
matrices can be found in the Expokit package [52[ • With the parameterization of Eq. (IA.1I) , the integral in Eq. (|A.4|) 
is elementary. 

A more accurate flavor evolution (beyond the constant-density approximation) can be obtained by applying the 
Magnus expansion [53], where the evolution operator is written as the exponential of an operator series, namely 



T(t) = exp 



(A.6) 



with 



fti(t) = -i [ dtiUx, 
Jo 

n 2 (t) = -\ I dh / 1 dt 2 [Hi,H 2 ] , 
* Jo Jo 

fia(t) = \ I dh f 1 dt 2 [ 2 dt 3 ([Hi,[H 2 , Ha]] + [Ha, [%2,Ui]]) , (A.7) 
Jo Jo Jo 



and so on, where we have used the shorthand Hi = H(ti). At first order, the Magnus expansion returns Eq. (IA.3I) . 
At second order, it is \K\,V. 2 \ = [Hy, V(x2) — V(ii)]. Integrating by part, one obtains 

T (Pk-iPk) = exp [-iHf ■ D k ] , (A.8) 

with 

Hf =U v + V k + i [H v , M k ] , (A.9) 

where 

M k = r dx V(x) (x - Xk -\ +Xk ) (A.10) 



D 



Xk-l 



is the "first moment" of the matter potential around the trajectory midpoint inside the fc-th shell. By using Eq. (|A.5|) 
and the parameterization in Eq. (IA.1I) , the integral in Eq. (|A.10|) is elementary. 

Concerning the flavor evolution of atmospheric neutrinos, we have adopted the second-order Magnus expansion for 
generic (real or complex) Hamiltonian, and we have checked that this approximation retains all the advantages of 
a fast analytical solution, without introducing significant differences with respect to the more accurate (but slower) 
numerical integration along the Earth density profile. We have also checked that our codes reproduce well the 
oscillograms discussed in [54[ (not shown). 
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